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We analyse a system in which, due to entanglement between the spin and spatial degrees of 
freedom, the reduced transmitted state has the shape of the freely propagating pulse translated in 
the complex co-ordinate plane. In the case an apparently "superluminal" advancement of the pulse 
the delay amplitude distribution is found to be a peculiar approximation to the Dirac 5-function, 
and the transmission coefficient exhibits a well-defined super-oscillatory window. Analogies with 
potential tunnelling and the Wheeler's delayed choice experiment are highlighted. 
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Introduction. Quantum speed up effect known since 
early thirties of last century allows the maximum of 
a wavepacket, transmitted across a classically forbid- 
den region, arrive at a detector earlier than that of a 
freely propagating one. It has been predicted and ob- 
served for various systems such as potential barriers, 
semi-transparent mirrors, refraction of light, microwaves 
in undersized wave guides and fast-light materials [1] , [2] 
A transmission appears "superluminal" if one uses the 
advancement of the transmitted peak to predict the time 
spent by the particle, photon or electromagnetic pulse 
inside a barrier, scatterer or a waveguide. It is generally 
understood that such a prediction is incorrect since the 
initial pulse, usually greatly reduced, undergoes in the 
process a severe reshaping. 

Superoscillatory functions which, in a limited region, 
oscillate with a frequency outside the support of their 
Fourier spectrum were introduced in [5] and extensively 
studied in [3] . These functions have been recently applied 
for improving optical resolution beyond the diffraction 
limit [5] . Authors of Ref. [5] have established a connec- 
tion between superoscillatory behaviour and anomalous 
tunnelling times (weak values) which occur in "super- 
luminal" signal transmission. The analogy between ap- 
parent "superluminality" and weak measurements was 
further explored in Ref. [7], and additional work on the 
weak values and superoscillations can be found in [5] . 

The purpose of this paper is to propose and analyse 
a system in which a similar speed up effect occurs due 
to entanglement between the spin and spatial degrees 
of freedom in a magnetic field so that the transmitted 
state is essentially a reduced copy of the freely propa- 
gating one, with an additional coordinate shift. We will 
show that in the case of a significant advancement the 
delay amplitude distribution [7J (DAD) is a wildly oscil- 
latory function with a finite support f2 which approxi- 
mates a Dirac (^-function with support outside £1. In the 
momentum space, the transmission coefficient exhibits 
a superoscillatory behaviour within a well defined band. 



We will also demonstrate that manipulations of the par- 
ticle's spin states allows to choose between the regimes 
where the particle spend a known amount of time in the 
magnetic field and where, in the spirit of wave-particle 
duality, such a duration cannot be defined. 

Quantum speed up effect via entanglement. We con- 
sider a particle of a unit mass /i = 1 equipped with a 
2K + 1-component magnetic moment (spin) . The parti- 
cle is described by a wave packet with a mean momentum 
Po, which propagates without spreading in the time inter- 
val of interest. Thus, its freely propagating wave function 
is given by 

(x\V(t)) = exp(ip Q x - iplt/2)G(x - po*) | a) (1) 
where G is the envelope and \a) = 
Em=-K a m \m)/y/N(a) 

K 

N(a)= M 2 ( 2 ) 

m=-K 

is the initial state of the spin written here in terms of 
its components, m, along the z-axis. Next we make the 
particle pass through a constant magnetic field in the z- 
direction created in the region < x < d so that the 
m-th component of the spin wave function encounters 
there an additional rectangular potential of the magni- 
tude mwi (lol is the Larmor frequency), a well if m < 
and a barrier for m > 0. If the particle is sufficiently fast, 
Pq/2 >> Kujl, the reflection off the field's edges can be 
neglected, and upon traversing the field the m-th compo- 
nent of the wave function will be advanced or lagging be- 
hind the freely propagating pulse by mAx, Ax = oj^cI/pq. 
On exit from the magnetic field the wave function con- 
tains a superposition of shifted wave packet's envelopes, 
with shifts ranging from — KAx to KAx, 

K 

(x\^(t)) = exp(ip x - iplt/2) a m x (3) 

7TV——K 

exp(—imui J d/po)G(x — pot — mAx)\m) j\f N(a). 
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For simplicity we will choose an initial spin state with 
a m — for all m > 0, so that none of the advanced 
shapes enter the sum ([3| , and post-select the particle spin 
in some known final state \b) = J2m=-K b m \m) / \f N (b) . 
We assume that such a post-selection can be performed 
as measuring the projector on \b), P(b), e.g., by making 
the particle pass through an additional polariser. On 
exit from the polariser the envelope of the particle's state 
(b|*(t)) becomes (X = x - p Q t): 



G(X) 



where 



G(X - x')r)(x')dx'/^N(a)N(b) 



K 



i](x) = ^ VmS(x + mAx), 
ra=0 

Vm = exp(-imw L d/p Q )a m b* 



(4) 



(5) 



We wish the envelope of the transmitted pulse to be an 
accurate (although possibly reduced) copy of the free en- 
velope translated in space by a distance a, 



G(x, t) m G(X - a)/y/N(a)N(b). 



(6) 



Choosing a > gives the impression that the particle has 
been sped up while passing through the set up containing 
the the magnetic field. 

Coordinate space: quasi-Dirac distributions. Equation 
^ requires the delay amplitude distribution (DAD) rj(x) 
in Eq.Q to be Dirac's S(X — a), something that cannot 
be achieved with the proposed set up. We can, however, 
use the freedom in choosing the spin's initial and final 
states | a) and \b) to insure that the normalisation and 
the first K moments of r](x) equal those of 6(X — a), 



x n rj(x)dx 



0, 1, ...K. 



(7) 



This is equivalent to K + 1 linear equations for the un- 
knowns rj m , 



-K 



Antrim = a", n = 0, 1, ...K 



(8) 



with a van der Monde matrix |9 A n , m = (mAx) n . 
Note that for i] m > Eq. ([7]) has a unique solution 
i](x) = 8(x — a). With no such restriction, other non- 
trivial solutions of Eqs.(|8]) are possible. Solving Eqs. ^ 

in a usual manner [10] yields (Y[ J= o indicates the prod- 
uct over all j ^ m) 



rj m (a/Ax) = (-1) 



m\(K — m)\ 



(9) 



With the quantities rj m (a) defined in Eq. Q we have a 
somewhat unusual mathematical object, an alternating 



distribution rj(x) with a finite support inside [— KAx, 0], 
which acts as the Dirac's 8(x — a) with support at an 
arbitrary a, positive or negative, in the space of all poly- 
nomials of order no higher than K. We may also expect 
it to have the same effect on any function G(x) in a do- 
main where its Taylor series can be truncated after K 
terms, 

j G(x)ri(x)dx w G (n) /n\ f x n r 1 (x)dx = (10) 

"* n=0 J 

K 

^G (n) a7n! « G(a), 



n=0 



although to check the validity of Eq. ( 10 1 we need, in 



principle, also analyse the behaviour of x n with n > K 
(FigjTJ). We will refer to an rj(x) with the above prop- 




FIG. 1: The moments of the DAD Q with K = 1, 15 and 30 
for a = a/(KAx) = 4. 

erties as a quasi-Dirac distribution (of order K). One 
example of such a distribution (of order 2) is the re- 
stricted path integral for the quantum traversal time [11] , 
which prompted Baz' [12] who compared its first and sec- 
ond moments, to conclude (erroneously |llj ) that the 
duration of an elastic collision time has a unique well 
defined value. We note further that a delay belonging 
to the discrete "spectrum" of available shifts, —MAx, 
— K < A I < 0, is produced by choosing either a m — 5 m M 
or b rn — S m M, and Eq. ^ yields, as it should, r] m — S m M 
(see Fig|2^). Choosing a delay between two available 
shifts, say, [M + 1/2) Ax, produces an alternating dis- 
tribution rj m plotted in Fig(2]3. Finally, an attempt to 
achieve a significant speed up (delay) outside the support 
of rj{x), \a/KAx\ > 1, yields r\ m which have very large 
absolute values, yet sum to unity as required by the first 
of Eqs. ([8]). The values r\ m for a/KAx = 4 are shown 
in Fig(2|:, where we had to use David Bailey's multiple 
precision program |13j to compute them to a sufficiently 
high accuracy. 

Probability of success and limits on the accuracy of 
transmission. Extremely large values of r\ m in Fig|2j; 
suggest that the proportion of the particles experiencing 
a significant advancement or delay will be very small as 
post-selection in the final state \b) will be unlikely to suc- 
ceed. Indeed, provided the envelope of the original pulse 
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FIG. 2: The DAD ^ for K = 30 and a/ Ax = -15, 
and 120 



-15.5 




10 -5 5 1 

X/(KAX) 



0-4 -2 2 4 

(P-Po)(Kax) 



is reproduced without distortion (i.e., Eq. ^ holds), 
the norm of the transmitted state is given by P(z) rs 
l/N(a)N(b) = [T, m= -K^T,° m =-K\Vm\ 2 /z m }-\ z m = 
\a m \ 2 where we have used \b m \ = \r] m \/\a m \, implied 
in Eq. ([5]). Maximazing P(z) with respect to z m , 
m = -K, ...0 yields \a m \ 2 = C\r/ m \, \b m \ 2 = \rj m \/C, 
where C is a constant independent of m so that the best 
probability of successful post-selection is given by 



FIG. 3: a) The transmitted envelope (multiplied by the fac- 
tor l/VP beat in Eq.|ll])) for K = 30, a/KAx = 2 and 



a = a/(KAx) = 4. Also shown by the dashed line is G(X — a) 
in the r.h.s. of Eq.|6|; b) same as a) except for a — 4.5; c) 
same as a) except for a = 3.5 + 2i; d) Re(T(p)) vs. p (solid) 
for the case shown in a). Also shown are sin(-ap) (dashed) 
and A(p) normalised to a unit height (thick solid). Vertical 
lines and arrows indicate the boundaries of the supersocilla- 
tory window; e) same as d) but for the case shown in b); f) 
\T(p)/exp(~ iaip + ct2p)\ (solid) and A(p) (thick solid) for the 
case shown in c). 



P best (a) = l/( ]T 



(11) amplitude T(p) is the Fourier transform of the DAD r)(x), 



-K 



As it has been shown above, for \a/K Ax \ > 1, 
Y^n^Khnl is, unlike Yfm^-K 1 !™ = 1, a very large 
number, thus making P best (a) << 1. 
For given states \a) and the success of advancing or 
delaying the free pulse depends on the shift a and the 
shape of the initial envelope G(X), which we choose to 
be a Gaussian with a coordinate width a. 



G(X) = (2/ 7 r C r 2 ) 1 / 4 exp(-A: 2 /o- 2 )- 



(12) 



The shapes of the transmitted pulse are shown in Figs(3] 
a,b and c. It is seen that an increase of the shift a results 
in the rapid onset of oscillations behind the peak of the 
transmitted pulse. 

Momentum space: superoscillations, real and complex- 
valued. To examine relation between the quasi-Dirac dis- 
tributions of Sect. 3 and superoscillatory functions intro- 
duced in jS] we expand the envelope G{X) in plane waves, 



G(X)= / T(p)A(p)exp(ipX)d P /^N(a)N(b), (13) 



where G{x) = J A(p) exp(ipx)dp and the transmission 



T(P) 



K 

I r](x) exp(—ipx)dx = rj m {a) exp(impAx) 
J m=0 



(14) 

We note that for r/(x) = 6(x — a), Eq. (14 1 would yield 
T(p) = exp(-ipa). With rj{x) being only a K-th order 
approximation to the Dirac 5(x — a) we can expect such 
a behaviour in a limited range of p's only. For an ad- 
vancement a > 0, T(p) which builds up from exponen- 
tial with non-negative frequencies mAx, m = 0,1. ..K, 
will behave there as an exponential with a negative fre- 
quency —a, i.e., will be superoscillatory. This superoscil- 
latory window (band) is essentially the region around 
p = where both T(p) and exp(— iap) are correctly 
reproduced by the first K terms of their Taylor series. 
Writing T{p) = J2n=o d^T(0)p n /n\, and recalling from 
Eqs.([7]) and Q 



6£T(0) = (■ 



-i) k x k 



k„k 



H) fc a 



< k < K 



(15) 



we find that the A"-th term is negligible provided 



(|p||a|)*/Jin « 1, 



for, 

\P\ < K/(e\a\), 



(16) 



where we have used the Stirling formula for the factorial. 
We can also choose the shift a to be complex valued, 
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a = ai+ia2, in which case one with a complex frequency, 
we have T(p) w exp(— ipai + pa2), and find the free en- 
velope in Eq.([6| analytically continued into the complex 
coordinate plane. In general, in order to avoid distor- 
tion we need the momentum distribution of the initial 
wavepacket A(p) to fit into the supcroscillatory band of 
the transmission coefficient, beyond which T(j>) exhibits 
rapid (polynomial) growth similar to that reported in 
Refs. |Sj,|S] (Fig 3L For a Gaussian envelope ( |l2| ) this 
requires (c.f. Eq.(16)) 2/<r << K/(e\a\) which puts a re- 



striction on the minimal width of the wave packet which 
can be translated without distortion. 

Conclusions and discussion. In summary, we have 
proposed a system for which the transmitted state is a 
reduced copy of the freely propagating one, translated 
along an arbitrary direction into the complex coordinate 
plane. The shift is achieved by splitting the initial pulse 
into several copies delayed relative to free propagation. 
The system can be described in terms of a transmission 
amplitude T(p) or, alternatively, by its Fourier trans- 
form, the delay amplitude distribution (DAD) i](x). The 
DAD contains information about the spectrum of virtual 
delays experienced by the particle in the magnetic field 
and, therefore, about causality. For a real translation the 
DAD is a finite order approximation to a Dirac <5-function 
and, for a complex valued shift, an analytic continuation 
of such an approximation. Accordingly, T(j>) exhibits a 
superoscillatory band whose width limits the size of the 
shift and the width of the initial pulse. The bandwidth 
depends on the number of first moments of the DAD 
(momentum derivatives of T{p) at p = 0) which satisfy 
x 11 Ria". A similar analysis can be applied to any system 
described by a transmission amplitude which is analytic 
in the complex momentum plane, e.g., tunnelling across a 
potential barrier [7] , where the shift is determined by the 
barrier width and the DAD, like that in Eq. Q contains 
only delayed contributions [7]. 

It is worth noting again that certain contradictions 
in the analysis of the so-called tunnelling time problem 
[1] arise in an attempt to deduce the time r a particle 
spends in the scatterer from its delay AT = —a/po in 
arriving at a remote detector. This suggests that the 
time spent in the scatterer is r = d/pg — a/po and 
the notion of apparent "superluminality" arises if r is 
made to be less than d/c. Our model readily demon- 
strates that there is no simple relation between AT and 
t. Indeed, between the field and the polariser the prob- 
ability density corresponds, as it should, to K real de- 
lays, P{x) = Y? m =-K\a m \ 2 \G{x~p t-m&x)\ 2 . The 
advanced shape (JsJ) is formed upon passing a (possibly 
remote) polariser, by an interference mechanism analo- 
gous to that responsible for anomalous values obtained 
in weak measurements |14j . There is also a similarity to 
the Wheeler's delayed choice experiment [15] ■ Just like 
Wheeler's photon may behave as a particle or as a wave, 
depending on which observation is made after the event, 



the particle spends in the magnetic field a well defined 
duration T = [d — mAx) /po provided the polariser is set 
to select one of the states \m). If, on the other hand, 
the post-selection is performed in a linear combination 
of more than one of such states, this duration cannot be 
defined and the experiment probes the wave nature of 
the transmitted particle. 

Acknowledgements: DS and RSM are grateful to Max- 
Plank-Institute for Physics of Complex Systems (Dres- 
den) for hospitality and financial support; RMS acknowl- 
edges Ministerio de Education y Ciencia, Plan National 
for support under grant No. FIS2007-64018. 



[1] E. H. Hauge and J. A. Stoevneng, Rev. Mod. Phys. 61, 
917 (1989); C. A. A. de Carvalho, H. M. Nussenzweig, 
Rev. Mod. Phys. 364, 83 (2002); V. S. Olkhovsky, E. 
Recami and J. Jakiel, Rev. Mod. Phys. 398, 133 (2004). 

[2] Y.P.Wang and D.L. Zhang, Phys. Rev .A 52, 2597 

(1995) ; Y.Japha and G.Kurizki, Phys. Rev .A 53, 586 

(1996) ; X.Chen and C.F.Li, Eur.Phys.Lett, 82, 30009 
(2008). 

[3] M. V. Berry, J. Phys. A, 27, L391 (1994). 

[4] A. Kempf, J. Math. Phys., 41, 2360 (2000); A. Kempf 
and P.J.S Ferreira, J. Phys. A: Math. Gen., 37, 12067 
(2004); M. S. Calder and A. Kempf, J. Math. Phys., 46, 
012101 (2005); P. J. S .Ferreira, A. Kempf and M. J. C. 
S. Reis, J. Phys. A, 27, L391 (2007). 

[5] M. V. Berry and S. Popescu, J. Phys. A, 39, 6965 (2006); 
F. M. Huang, Y. Chen, F. J. Garcia de Abajo and N. I. 
Zheludev, J. Opt. A: Pure Appl. Opt. 9, S285 (2007); M. 
R. Dennis, A. C. Hamilton and J. Courtial, Opt. Lett., 
33, 2976 (2008); N. I. Zheludev, Nature (materials), 7, 
420 (2008); M. V. Berry and M. R. Dennis, J. Phys. A: 
Math. Theor. 42, 022003 (2009). 

[6] Y. Aharonov, N. Erez and B. Reznik, Phys. Rev. A, 65, 
052124 (2002). 

[7] D. Sokolovski, A. Z. Msezane and V. R. Shaginyan, Phys. 

Rev. A, 71, 064103 (2005). 
[8] J. Tollaksen, J. Phys.: Conf. Series, 70, 012016 (2007). 
[9] A. Ostaszewski, Advanced Mathematical Methods (Cam- 
bridge University Press, 1990). 

[10] Consider the determinant A m (a) of the matrix A in 
which the m-th column is replaced by (1, a, a 2 , ...a K ). 
A m (a), which is a polynomial of order K in a, must 
vanish for any a — —nAx, n 7^ m, and, therefore, can 
be written as A m (a) = Iln^mC + n Ax) x const. Also, 
we have detA = Aj(-jAx), for any < j < K , and the 
ratio A m (a)/detA yields Eq.|9]|. 

[11] D. Sokolovski, Phys. Rev. A, 76, 042125 (2007). 

[12] A. I. Baz': Yad. Fiz. 4, 182 (1967) [Sov.J.Nucl.Phys, 5, 
635 (1967)]. 

[13] D. H. Bailey, ACM Transactions on Mathematical Sot- 

fware, 19, 288 (1993). 
[14] Y. Aharonov, D. Albert and L. Vaidman, Phys. Rev. A, 

60, 1351 (1988). 
[15] J. A. Wheeler, in Mathematical Foundations of Quantum 

Theory, edited by A. R. Marlow, Academic Press, New 

York, 1978, p. 13. 



